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O ! ABSTRACT 
^ . 

I Aims. We investigate how ellipticity, asymmetries and substructures separately affect the ability of galaxy clusters to produce strong tensing 
. events, i.e. gravitational arcs, and how they influence the arc morphologies and fluxes. This is important for those studies aiming, for example, 

CIh' at constraining cosmological parameters from statistical tensing, or at determining the inner structure of galaxy clusters through gravitational 

^ arcs. 

» I Methods. We do so by creating two-dimensional smoothed, dilferently elliptical and asymmetric versions of some numerical models. By 
' subtracting these smoothed mass distributions from the corresponding numerical maps and by gradually smoothing the residuals before re- 
^ adding them to the clusters, we are able to see how the lensing properties of the clusters react to even small modification of the cluster 
|L* morphology. We study in particular by how much ellipticity, asymmetries and substructures contribute to the strong lensing cross sections of 
. clusters. We also investigate how cluster substructures affect the morphological properties of gravitational arcs, their positions and fluxes. 

■ Results. On average, we find that the contributions of ellipticity, asymmetries and substructures amount to ~ 40%, ~ 10% and ~ 30% of the 
$_( ' total strong lensing cross section, respectively. However, our analysis shows that substructures play a more important role in less elliptical and 
asymmetric clusters, even if located at large distances from the cluster centres (~ Ih 'Mpc). Conversely, their elfect is less important in highly 
asymmetric lenses. The morphology, position and flux of individual arcs are strongly alfected by the presence of substructures in the clusters. 
Removing substructures on spatial scales < 50/j"'kpc, roughly corresponding to mass scales < 5 x 10'"/;"' Mq, alters the image multiplicity 
of ~ 35% of the sources used in the simulations and causes position shifts larger than 5" for ~ 40% of the arcs longer than 5". 
Conclusions. We conclude that any model for cluster lens cannot neglect the effects of ellipticity, asymmetries and substructures. On the 
other hand, the high sensitivity of gravitational arcs to deviations from regular, smooth and symmetric mass distributions suggests that strong 
gravitational lensing is potentially a powerful tool to measure the level of substructures and asymmetries in clusters. 

Key words, gravitational lensing - galaxies: clusters, dark matter 



1. Introduction 

Thanks to the improvements in the quality and in the depth 
of astronomical observations, in particular from space, an in- 
creasing number of gravitational arcs has recently been dis- 
covered near the cent res of many galaxy clusters (see e.g. 
iBroadhurst et alJ l20051. Since the appearance of these images 
reflects the shape of the gravitational potential which is respon- 
sible for their large distortions, strong lensing is, in principle, a 
very powerful tool for investigating how the matter, in partic- 
ular the dark component, is distributed in the inner regions of 
cluster lenses. 



Determining the inner structure of galaxy clusters is one 
of the major goals in cosmology, because it should allow us 
to set important constraints on the growth of the cosmic struc- 
tures in the Universe. Moreover, constraining the mass distribu- 
tion in the centre of dark matter halos has become increasingly 
important in the recent years, since observations of the dy- 
namics of stars in galaxy-sized systems revealed the presence 
of a potential problem within the Cold-Dark-Matter (CDM) 
scenario. While numerical simulations in this cosmological 
framework predict that dark matter halos in a large range of 
masses should develop density profiles characterised by an in- 
ner cusp, observations of the rotation curves of dwarf and low- 
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surface-brightness galax ies sugge st tha t these objects rather 
have flat density profile s (Flores & Primack 1994; Moore 1994; 
Burkert* '1995'; 'fiurkert & Silk 



1997; McGaugh & de Blok 
1998; Dalcan ton & Bernstein 2000; Firmani et al. 2001). 



While the centres of galaxies are dominated by stars, which 
renders it extremely complicated to derive constraints on the 
distribution of their dark matter, galaxy clusters are an alterna- 
tive and, in many respects, preferable class of objects for test- 
ing the predictions of the CDM model. In fact, several authors 
already tried to investigate the inner structure of these large 
systems, using and often combining several kinds of observa- 
tions. Apart from lensing, the gravitational potential of galaxy 
clusters can be traced with several other methods, for example 
through the emission in the X-ray band by the hot intra-cluster 
gas. However, while gravitational lensing directly probes the 
matter content of these objects, the other techniques usually 
rely on some strong assumptions about their dynamical state 
and the interaction between their baryonic and dark matter For 
example, it must be often assumed that the gas is in hydrostatic 
equilibrium within the dark matter potential well and that the 
system is spherically symmetric. 

Some ambiguous results were found when comparing 
the constraints on the inner structure of clusters as ob- 
tained from X-ray and lensing observations. First, masses 
estimated from strong lensing are usually larger by a fac- 
tor o f 2-3 than the rnasses obtained from X-ray observa- 
tions jChen et 31.1120031 lOta et al.ll2004h . Deviations from ax- 
ial symmetry and substructures are known to be important 
factor s in strong lensing mass estimates (see e.g. |b artelmanij 
'l995; Bartelmann & Steinmetz 1996; lMeneghetti et al. 2003b; 
Qguri et al...2005.' .Gavazzi.2005 ). Second, the constraints on 



the inner slope of the density profiles se em to be compat - 
ible with a wide range of inner s lopes (Ettori et al. 2002; 



Lewis et al.' '2003', 'Arabadiis et al.' '2002'; Sa nd et alj 12003 : 
Bartelmann & Meneghetti 2004; Gavazzi 2005). 

Apart from the above-mentioned uncertainties affecting the 
X-ray measurements, strong lensing observations also have 
several potential weaknesses. First of all, arcs are relatively rare 
events. Frequently, all the constraints which can be set on the 
inner structure of clusters via strong lensing depend on a sin- 
gle or on a small number of arcs and arclets observed near the 
cluster core. Second, arcs are the result of highly non-linear 
efifects. This implies that their occurrence and their morpho- 
logical properties are very sensitive to ellipticity, asymmetries 
and substructures of the cluster matter distribution. 

Reversing the problem, this means that, in order to re- 
liably describe the strong lensing properties of galaxy clus- 
ters, all of these effects must be taken into account. Fitting 
the positions and the morphology of gravitational arcs to de- 
rive the underlying mass distributions of the lensing clusters, 
usually requires to build models with multiple mass compo- 
nents, each of whi ch is charac terised b y its ellipticity and 
1^ 



orientation (see e.g. 'K neib et aD>1993; 



iBro adhurst et al. 2005). Even describing the cluster lens popu- 
lation in a statistical way re quires to use realistic cluster mod- 
els tMeneghetti et alJEoOoi i2003b.a; .OgurL2002; .Qguri et alJ 
l2003UDalal et al.l2005HHennawi et alJ2005l) . 



Despite the fact that the importance of ellipticity, asym- 
metries and substructures for strong lensing appears clearly in 
many previous studies, many questions still remain. For exam- 
ple, what is the typical scale of substructures which contribute 
significantly to the strong lensing ability of a cluster? Where 
are they located within the clusters? What is the relative impor- 
tance of asymmetries compared to ellipticity? Moreover, how 
do substructures influence the appearance of giant arcs? All of 
these open problems are important for those studies aiming at 
constraining cosmological parameters from statistical lensing, 
or at determining the inner structure of galaxy clusters through 
gravitational arcs. 

This paper aims at answering to these questions. To do so, 
we quantify the impact of ellipticity, asymmetries and substruc- 
tures by creating differently smoothed models of the projected 
mass distributions of some numerical clusters. We gradually 
move from one smoothed model to another through a sequence 
of intermediate steps. 

The plan of the paper is as follows. In Sect.|2l we discuss 
the characteristics of the numerically simulated clusters that 
we use in this study; in Sect. |3j we explain how ray-tracing 
simulations are carried out; Sect. |4] illustrates how we obtain 
smoothed versions of the numerical clusters; in Sect.|5l we sug- 
gest a method to quantify the amount of substructures, asym- 
metry and ellipticity of the cluster lenses, based on multipole 
expansions of their surface density fields; Sect. |6lis dedicated 
to the discussion of the results of our analysis. Finally, we sum- 
marise our conclusions in Sect0 

2. Numerical models 

The cluster sample used in this paper is made of five massive 
dark matter halos. One of them, labelled g8hi, was simulated 
with very high mass resolution, but contains only dark-matter. 
The others, the clusters gl, §8, ^51 and g72 have lower mass 
resolution but are obtained from hydro-dynamical simulations 
which also include gas. 

The halos we use here are massive objects with masses 
8.1 xW^h-^MQ iglT), 8.6x lO'^/i-iMo {g5\), 1 AxW^h-^M^ 
(g\) and 1.8 x 10^^ h-^M^ (g8 and ^8hr) at z = 0.3. We have 
chosen this redshift because it is close to where the strong lens- 
ing efficien cy of clusters is the largest for sources at Zs ^ 1 
(iLi etalj2 005). 

The clusters were extracted from a cosmological simulation 
with a box-size of 479/1 'Mpc of a flat ACDM model with 
Qq ^ 0.3, h = 0.7, 0-8 = 0.9, and Ob = 0.04 (see Yoshida e tall 
2001) . Using the "Z oomed Initial Conditions" (ZIC) technique 
( Tormen et alJll997i) . they were re-simulated with higher mass 
and force resolution by populating their Lagrangian volumes 
in the initial domain with more particles, appropriately adding 
small-scale power. The initial displacements are generated us- 
ing a "glass" distribution (White 1996) for the Lagrangian par- 
ticles. The re-sim ulations wer e carried out with the Tree-SPH 
code GADGET-2 dSpringel et al..200k.SpringeL2005>) . For the 
low resolution clusters, the simulations started with a gravi- 
tational softening length fixed at e = 30.0/; 'kpc comoving 
(Plummer-equivalent) and switch to a physical softening length 
of e = 5.0/j^'kpc at 1 H-z = 6. 
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The particle masses are moM = 1.13 x 10^ /j^'Mq and 
'Wgas = 1.7 X 10^ /z~'Mo for the dark matter and gas parti- 
cles, respectively. For the high-resolution cluster ^Sh, the parti- 
cle mass is ;«dm = 2.0 x 10^/i"'Mq and the softening was set to 
half of the value used for the low resolution runs. Its virial re- 
gion at z = 0.3 contains more than nine million particles, which 
allow us to well resolve substructures on scales down to those 
of galaxies. Despite the lower mass resolution with respect to 
gShi, the other low resolution clusters also contain several mil- 
lion particles within their virial radii. 

To introduce gas into the high-resolution regions of the low- 
resolution clusters, each particle in a control run containing 
only dark matter was split into a gas and a dark matter particle. 
These were displaced by half the original mean inter-particle 
distance, so that the centre-of-mass and the momentum were 
conserved. 

Selection of the initial region was done with an iterative 
process involving several low-resolution, dissipation-less re- 
simulations to optimise the simulated volume. The iterative 
cleaning process ensures that all these haloes are free of con- 
taminating boundary effects up to at least 3 to 5 times the virial 
radius. 

The simulations including gas particles follow only the non 
radiative evolution of the intra-cluster medium. More sophisti- 
cated versions of these clusters, where radiative cooling, heat- 
ing by a UV background, and a treatment of the star forma- 
tion and feedback processes were included exist and their lens- 
ing^pro perties have been studied in detail by Puchwein et al. 

The cluster gShr is in principle a higher-resolution, dark- 
matter only version of the cluster §8, which was simulated 
with non-radiative gas physics. Nevertheless the two objects 
can only be compared statistically. Indeed, the introduction of 
the gas component as well as the increment of the mass res- 
olution introduce small perturbations to the initial conditions, 
which lead to slightly different time evolutions of the simulated 
halos. Furthermore, also the presence of gas and its drag due 
to pressure lead to significant changes in the assembly of the 
halo. A detailed discussion of such differences can be found in 
[Puchwein et al. (2005). For this reason, the lensing properties 
of gS and ^8;,^ at z = 0.3 are not directly comparable. 

3. Lensing simulations 

Ray-tracing simulations are carried out using the technique de- 
scribed in detail in several e arlier papers (e.g. Bartelmann et al. 
ll998HMeneghetti et all 2000). 

We select a cube of 6/; 'Mpc comoving side length, cen- 
tred on the halo centre and containing the high-density region 
of the cluster The particles in this cube are used for produc- 
ing a three-dimensional density field, by interpolating their po- 
sition on a gri d of 1024-^ cells using the T riangular Shaped 
Cloud method ( iHocknev & Eastwoo Then, we project 

the three-dimensional density field along the coordinate axes, 
obtaining three surface density maps used as lens planes in 
the following lensing simulations. 

The lensing simulations are performed by tracing a bundle 
of 2048 X 2048 light rays through a regular grid, covering the 
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central sixteenth of the lens plane. This choice is driven by the 
necessity to study in detail the central region of the clusters, 
where critical curves form, taking into account the contribution 
from the surrounding mass distribution to the deflection angle 
of each ray. 

Deflection angles on the ray grid are computed using the 
method described in Meneghetti et al. (2000). We first define a 
grid of 256 x 256 "test" rays, for each of which the deflection 
angle is calculated by directly summing the contributions from 
all cells on the surface density map 

where A is the area of one pixel on the surface density map 
and Xh,k and Xij are the positions on the lens plane of the "test" 
ray (h, k) and of th e surface density element (/, j). Following 
IWambsganss et all (|l998), we avoid the divergence when the 
distance between a light ray and the density grid-point is zero 
by shifting the "test" ray grid by half-cells in both direc- 
tions with respect to the grid on which the surface density is 
given. We then determine the deflection angle of each of the 
2048 X 2048 light rays by bi-cubic interpolation between the 
four nearest test rays. 

The position y of each ray on the source plane is calculated 
by applying the lens equation. If y and x are the angular po- 
sitions of source and image from an arbitrarily defined optical 
axis passing through the observer and perpendicular to the lens 
and source planes, this is written as 

y = x- j^aix) , (2) 

where D\s and are the angular diameter distances between 
the lens and the source planes and between the observer and 
the source plane, respectively. 

Then, a large number of sources is distributed on the source 
plane. We place this plane at redshift Zs = 2. Keeping all 
sources at the same redshift is an approximation justified for 
the purposes of the present case study, b ut the recent detections 
of ar cs in high-reds hift clusters ( Zaritskv & Gonzale j l2003l 
Glad ders et al]l2003l) indicate that more realistic simulations 
will have to account for a wide source redshift distribution. 

The sources are elliptical with axis ratios randomly drawn 
from [0.5, 1]. Their equivalent diameter (the diameter of the 
circle enclosing the same area of the source) is re - 1". They 
are distributed on a region on the source plane corresponding 
to one quarter of the field of view where rays are traced. As 
in our earlier studies, we adopt an adaptive refinement tech- 
nique when placing sources on their plane. We first start with 
a coarse distribution of 32 x 32 sources and then increase the 
source number density towards the high-magnification regions 
of the source plane by adding sources on sub-grids whose reso- 
lution is increased towards the lens caustics. This increases the 
probability of producing long arcs and thus the numerical effi- 
ciency of the method. In order to compensate for this artificial 
source-density enhancement, we assign a statistical weight to 
each image for the following statistical analysis which is pro- 
portional to the area of the sub-grid cell on which the source 
was placed. 
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By collecting rays whose positions on the source plane 
lie within any single source, we reconstruct the images of 
background galaxies and measure their length and width. 
Our technique for i mage detection a nd classification was de- 
scribed in detail by Bartelmann & Weiss ( 199 41) and use d by 
Meneehetti etal.'(2000, 2001, 2003a ta): iiTorri e'tal]ll2004l ) and 
Meneehetti et al. (2005). The modifications recently suggested 
by IPuchwein et alJ ll2005ll for increasing the accuracy of the 
measurements of the arc properties have been included in our 
code. The simulation process ends in a catalogue of images 
which is subsequently analysed. 



4. Smoothed representations of the cluster 

We aim here at separating the effects of substructures, ellip- 
ticity and asymmetries on the strong lensing efficiency of our 
numerical clusters. We work in two dimensions, starting from 
the projected two-dimensional mass distribution of each lens. 
Since the lens' surface density profile plays a crucial role for 
strong lensing, we keep it fixed and vary only the shape of the 
density contours. 

As a first step, we construct a fiducial model for the smooth 
mass distribution of the lens. This is done by measuring the 
ellipticity and the position angle of the surface density con- 
tours as function of radius in the projected mass map. The pro- 
jected density is measured in circles of increasing radii x. We 
determine the quadrupole moments of the density distribution 
in each aperture, 

J d-^xl,(x) 

where Xc denotes the position of the cluster centre and the 
integrals are extended to the area enclosed by the aperture. Both 
the ellipticity e and the position angle <p of the iso-density con- 
tour corresponding to the chosen aperture radius are derived 
from the elements of the tensor 5, ;(r): 



e(x) = 



(5 11 +522)2 



4>(x) = -arctan- — 

Z i 11 - i 22 



(4) 
(5) 



The ellipticity and position angle profiles are used in com- 
bination with the mean surface density profile S(x) of the lens 
for constructing a smoothed surface density map, 

= nx,,^) , (6) 

where the equivalent radius x^ ^ is given by 

^ [xi cos <p(x) + X2 sin (p(x)]^ 



[l-e{x)] 



[-x\ sin (pix) + X2 cos (;A(x)]2[l + e(x)]) 



1/2 



(7) 



and e = 1 - Vl -e^lil + e) = 1 - VT^/ VT+l ^ e (for 
e<K 1). 



The resulting map conserves the mean surface density pro- 
file of the cluster and reproduces well the twist of its iso- 
density contours, i.e. the asymmetries of the projected mass 
distribution. This is shown in Fig. [2 The panel on the left col- 
umn shows the original surface density map for one projec- 
tion of the high-resolution cluster ^8hi . Contour levels start at ~ 
3.6 X 10'^ h Mq Mpc"2 and are spaced by 5 x 10'^ h Mq Mpc"^. 
The top-middle panel shows the smoothed map obtained from 
Eq. (|6j. In the rest of the paper, we will call this smoothed 
model the "asymmetric" model. The same colour scale and 
spacing of the contour levels as in the first panel are used. 
In the smoothed map, the substructures on all scales are re- 
moved and redistributed in elliptical shells around the cluster. 
Comparing the strong lensing properties of the original and of 
the smoothed map, we can therefore quantify the net effect of 
substructures on the cluster strong lensing efficiency. By sub- 
tracting the smoothed from the original map, we obtain a resid- 
ual map showing which substructures will not contribute to 
lensing after smoothing. We plot this residual map in the right 
panel in the upper row of Fig. [2 

Similarly, the eff'ects of other cluster properties can be sepa- 
rated. For example, we can remove asymmetries and deviations 
from a purely elliptical projected mass distribution by disabling 
the twist of the iso-density contours in our smoothing proce- 
dure. For doing so, we fix the ellipticity and the position angle 
to a constant value, e - ecrit and <p = 0ciit- We choose fciit and 
ipait to be those measured in the smallest aperture containing 
the cluster critical curves. A smoothed map of the cluster is cre- 
ated as explained earlier The results are shown in the middle 
panels of Fig. [2 A comparison of the lensing properties of this 
new "elliptical" model with those of the previously smoothed 
map allows us to quantify the effect of asymmetries, which 
are large-scale deviations from elliptical two-dimensional mass 
distributions. 

Finally, even the cluster ellipticity can be removed, still pre- 
serving the same mean density profile. This is easily done by 
inserting e = in Eq. (|6}. The resulting smoothed map and 
the residuals obtained by subtracting it from the original pro- 
jected mass map are shown in the bottom panels of FigQ] If 
we compare the lensing efficiency of such an "axially symmet- 
ric" model to that of the previously defined elliptical model, 
we quantify the effect of ellipticity on the cluster strong lensing 
properties. 

For each smoothing method, we simulate lensing of back- 
ground galaxies not only for the extreme cases of the totally 
smoothed maps but also for partially smoothed mass distribu- 
tions. Adding the residuals R to the smoothed map, the original 
surface density map of the cluster is indeed re-obtained. 



i:{x) = 5:™(a:) + R{x) . 



(8) 



Substructures of different sizes can be gradually removed from 
the cluster mass distribution by filtering the residual map with 
some filter function of varying width before re-adding it to the 
totally smoothed map. This can be done, for example, by con- 
volving the residual map R with a Gaussian G of width cr: 



R(x, cr) = R(x) * G{x, 0-) , 



(9) 
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Fig. 1. Different smoothing sequences of the projection along the x-axis of cluster §8hr- The original map is shown on the left. 
Top panels: smoothing with variable ellipticity and orientation of the iso-density contours; middle panels: smoothing assuming 
fixed ellipticity and position angle (see text for more details); bottom panels: smoothing assuming axial symmetry. From left to 
right, each row of panels show the smoothed map and the residuals, obtained by subtracting the smoothed from the original map. 
The contour levels in the density maps start at ~ 3.6 x 10'^ /z Mq Mpc - and are spaced by 5 x 10'^ /i Mq Mpc"^. The horizontal 
side length of each panel is 6 /z^'Mpc. 



G{x,cr) - 



(10) 



The width cr defines the scale of the substructures which will be 
filtered out of the mass map. Surface density maps with inter- 
mediate levels of substructures are finally obtained by adding 
the filtered residuals to the smoothed map, 



S(jc, 0-) = 2™(;t:) + R{x, cf) . 



(11) 



This procedure allows us to investigate what is the charac- 
teristic scale of substructures which are important for strong 
tensing. Moreover, comparing how the lensing properties of 
different clusters react to smoothing, we can quantify the im- 
pact of substructures in halos with different degrees of asym- 
metry and ellipticity. A sequence of smoothed versions of the 
same cluster model is shown in Fig|2 The smoothing length cr 
is 0,47, 141 and 470/? ' comoving kpc from the top left to the 
bottom right panel, respectively. 

Another important issue is to understand where the sub- 
structures must be located in order to have a significant im- 
pact on the strong lensing properties of clusters. To address 
this problem, we remove from the clusters the substructures 
located outside apertures of decreasing equivalent radius x^^^. 
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Fig. 2. Surface density maps of the same cluster projection as 
shown in Fig. ^ smoothed with increasing smoothing length 
from the top left to the bottom right panels. The background 
smoothed model is the one shown in the upper right panel of 
Fig-E The smoothing lengths in the four panels are 0,47, 141 
and 470/;" 'kpc comoving, respectively. The horizontal side 
length of each panel is 6 /; 'Mpc comoving. 
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Fig. 3. Surface density maps of the same cluster projection 
shown in Fig.^ suppressing the substructures in shells of de- 
creasing equivalent radius from the top left to the bottom right 
panels. The background smoothed model is the one shown in 
the upper right panel of Fig. The equivalent radii beyond 
which the substructures are removed are 1174,704,352 and 
235 /i 'kpc comoving, respectively. The horizontal side length 
of each panel is 6 /z 'Mpc comoving. 

Again, this is done by modifying the residuals of the smoothed 
maps. We multiply the residual map with the function. 



1 



exp 



( ^o-i] 



< I) 
(x,.4, > I) 



(12) 



where (Xcut = 100 /i 'kpc comoving and / is the equivalent 
distance beyond which the substructures are suppressed. The 
Gaussian tail of the window function was applied to avoid sud- 
den discontinuities in the surface density maps. 

Results of the removal of cluster substructures at different 
radii are shown in Fig.|3 From the top left to the bottom right 
panel, the cut-off scales I are 1 174, 704, 352 and 235 /i^'kpc co- 
moving, respectively. The residual maps filtered with the win- 
dow function (I12> were re-added to the totally smoothed asym- 
metric maps for obtaining surface density distributions with the 
desired level of substructures within a given equivalent radius. 

5. Quantifying the amount of substructures, 
asymmetry and ellipticity 

The variations of ellipticity and position angle of the iso- 
density contours are given by the functions e{x) and 
which were defined in the previous section. These are shown 
for the three projections of cluster ^8hr in Fig. |4] They illus- 
trate that the projection along the x-axis of this cluster, shown 
in Fig. n is the most elliptical at the relevant radii, with an 
ellipticity which grows from e ~ 0.4 to ~ 0.58 within the in- 
ner ~ 1 /j"'Mpc. The iso-density contours have almost constant 
orientation in this projection. The projection along the y-axis 
exhibits the largest variations of ellipticity in the central re- 
gion of the cluster, with e growing from ~ 0.25 to ~ 0.52. It 




Distance frcm the center [Mpc/h] 

Fig. 4. Variations of ellipticity (top panel) and position angle 
(bottom panel) of the iso-density contours of the three projec- 
tions of cluster gShi as function of the distance from the cluster 
centre. The cluster has a virial radius of ~ 2.7/! 'Mpc. 

is also characterised by a large twist of the iso-density con- 
tours, whose orientations change by up to ~ 20 degrees. When 
projected along the z-axis, the cluster appears more circular 
and with fairly constant ellipticity ranging between ~ 0.22 and 
~ 0.32. The twist of the iso-density contours is moderate within 
the inner 1 /; 'Mpc. 

We now quantify the amount of substructure within the nu- 
merical clusters by means of multipole expa nsions of their sur- 
face density maps dMeneghetti et alJ2003bl) . 

Briefly, starting from the particle positions in the numeri- 
cal simulations, we compute the surface density 2 at discrete 
radii x„ and position angles 6k taken from [0, 1.5] /i ' Mpc and 
[0, 27r], respectively. For any x„, each discrete sample of data 
S(jic„, Oil) is expanded into a Fourier series in the position angle, 



E(x„,0,) = J]5,(x„)e-'"'' , 

1=0 

where the coefficients S /(x„) 

oo 

5,(x„) = ^5:(x„,0,)e'"'S 



(13) 



(14) 



can be computed using fast-Fourier techniques. 
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We define the power spectrum Pi(x„) of the multi- 
pole expansion / as Pi( x„) - \Si(x„)f-. As discussed by 
iMeneghetti et alJ (l2003 lj). the amount of substructure, asym- 
metry and elUpticity in the mass distributions of the numeri- 
cally simulated cluster at any distance x„ from the main clump 
can be quantified by defining an integrated power Pciic(-x«) as 
the sum of the power spectra over all multipoles, from which 
the monopole is subtracted in order to remove the axially sym- 
metric contribution. 



power in substructures vs distance 



Pcuc(x„) = ^ Pl(x„) - Pu(x„) ■ 



(15) 



This quantity measures the deviation from a circular distribu- 
tion of the surface mass density at a given distance x„ from the 
cluster centre. 

In a fully analogous way, we can quantify the deviation 
from a purely elliptical surface mass density by subtracting 
from Pcirc the quadrupole term P2{xn): 



PelliXn) - Pciic(x„) - PliXn) . 



(16) 



Separating the effects of asymmetries and substructures by 
means of contributions to the multipoles is not an easy task. 
The two terms are mixed together in the multipole expansion. 
In this paper, asymmetries are assumed to be large-scale devi- 
ations from a purely elliptical mass distribution, which result 
in variations of ellipticity and position angle of the iso-density 
contours as functions of radius. Their contribution to the az- 
imuthal multipole expansion is mostly contained in the low-Z 
tail. On the other hand, for small x„ even relatively small sub- 
structures subtend large angles with respect to the cluster cen- 
tre. The dipole term itself contains a contribution from sub- 
structures. Assuming that substructures are localised lumps of 
matter whose angular extent is < n/l and then contribute to all 
multipoles of order larger than Z ~ 2, we quantify the amount of 
substructures at radius x„ in the cluster by means of the quantity 



P^uhiXn) = Peu(x„) - Pi(x„) , 



(17) 



where the power corresponding to the dipole, Pi(x„) has been 
subtracted from Peiiix„). 

Fig|5lshows the radial profiles of Psub for the three projec- 
tions along the x, y and z-axes of cluster g^hi-. Peaks along the 
curves indicate the presence of substructures. The amplitude of 
the peaks is a growing function of the mass of the substruc- 
tures. Clearly, in the projection along the x-axis massive lumps 
of matter are located at distances of ~ 550 and ~ 700/? 'kpc 
from the cluster centre. Substructures are abundant also at radii 
of ~ 1 and ~ 1.25 /z 'Mpc. This is visible in the left panel 
of Fig[2 Instead, the most dominant substructures in the pro- 
jection along the y-axis are located outside the region of ra- 
dius 1 /i 'Mpc. Only a relatively small peak is observed at 
~ 650/j"'kpc from the centre. Finally, when projected along 
the z-axis, the cluster contains a large amount of substructures 
at the distance of ~ 8OO/1 'kpc from the centre. Other peaks 
are located at radii > 1 /i^'Mpc. 
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Fig. 5. Power in substructures as a function of distance from 
the cluster centre for the three projections along the x-, y- and 
z-axes of cluster §8hr. 



6. Results 

In this section, we describe the lensing properties of the nu- 
merical clusters in the sample we have analysed and quantify 
the impact of ellipticity, asymmetries and substructures on their 
ability to produce arcs. 



6.1. Magnification and caustics 

The lensing properties of the two-dimensional mass distribu- 
tions generated using the previously explained methods can be 
easily determined using the standard ray-tracing technique de- 
scribed in Sect. 13 

The ability of galaxy clusters to produce strong lensing 
events is expected to reflect both the presence of substruc- 
tures embedded into their halos and the degree of ellipticity 
and asymmetry of their mass distributions. Indeed, all of these 
factors contribute to increase th e shear field of th e clusters. 
This was shown for example by 'Torri et al.' i a nd late r 
confirmed by Meneghetti et al. ( 2005 ) and Fedeli et all ( l2005l) . 
who found that the passage of substructures through the cluster 
cores can enhance the lensing cross section for the formation of 
arcs with large length-to -width ratios by orders of magnitude. 
IMeneghetti et al .'( 2003b) show that elliptical models with real- 
istic density profiles produce a number of arcs larger by a fac- 
tor of ten compared to axially-symmetric l enses with the sam e 
mass. Analogous results were obtained by Oguri et al. (20Q2I), 
who compared the lensing efficiency of triax ial and spherical ly 
symmetric halos, and more recently by Hen nawi et al ] |l2005l) . 

By smoothing the two-dimensional mass distribution of the 
clusters, both the levels of substructures and asymmetries are 
decreased. Thus, we expect their ability to produce highly dis- 
torted arcs to be somewhat reduced. This expectation is sup- 
ported by the fact that the regions of the lens plane where the 
tangential-to-radial magnification ratio exceeds a given thresh- 
old shrink significantly when the smoothing is applied. This 
is shown in Fig. |6l The map of the tangential-to-radial mag- 
nification ratio of the original cluster (top right panel) is com- 
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asymm. model 




sm. length = h"' kpc 

sm. length — 47 h"^ kpc 

sm. length = 141 h"' kpc 

sm. length = 470 h"' kpc 





Fig. 8. Lens caustics of the mass distribution of the cluster projected along the .JC-axis as obtained after smoothing using several 
smoothing lengths and assuming different background models: asymmetric model (left panel), elliptical model (central panel) 
and axially symmetric model (right panel). The side length of each panel is 500 /z 'kpc comoving, corresponding to ~ 1'. 




Fig. 6. Maps of the tangential-to-radial magnification ratio for 
the same cluster projection showed in the previous figures. 
From the top left to the bottom right panel, we show the maps 
corresponding to the original cluster and to three smoothed ver- 
sions of it: using the asymmetric, the elliptical and the axially 
symmetric background models. 



pared to those obtained by smoothing its surface density map 
using the asymmetric (top left panel), the elliptical (bottom left 
panel) and the axially symmetric (bottom right panel) back- 
ground models. The cumulative distributions of the pixel val- 
ues in these maps are displayed in Fig0 The probability of 
having pixels where the tangential-to-radial magnification ra- 
tio exceeds the minimal value decreases at least by a factor of 
two, due to removal of substructures, asymmetries and elliptic- 
ity. This does not imply that the cross section for arcs with large 
length-to-width ratio decreases by the same amount, since the 
excess of pixels with large tangential-to-radial magnification 
ratio in the unsmoothed map is in part due to isolated lumps of 
matter whose angular scale is similar or smaller than the angu- 
lar scale of the sources. 



original ciuster - 
asymm. mode 
eilipt. modei ■ 
ax. symm. modei 



Fig. 7. Cumulative distributions of tangential-to-radial magni- 
fications in the maps showed in Fig|6l 



By definition, the lensing cross section, which measures a 
cluster's ability to produce arcs, is an area encompassing the 
lens' caustics. Thus, the more extended the caustics are, the 
larger is generally the lensing cross section. In Fig.|8]we show 
how the caustic shape changes as the smoothing length is var- 
ied. Results are shown for each of the smoothing schemes ap- 
plied. As expected, the caustics shrink as the smoothing length 
increases. Comparable trends are found for the asymmetric and 
elliptical background models, for which the change of the caus- 
tic length is not dramatic. On the other hand, if the cluster 
surface density is gradually smoothed towards an axially sym- 
metric distribution, the evolution of the lens' caustics is much 
stronger 

Similar reductions of the caustic sizes are found when sup- 
pressing the substructures outside a given radius. This is shown 
in Fig. O Clearly, substructures at distances of the order of 
l/i 'Mpc already play a significant role for strong lensing. 
Although they are located far away from the cluster critical re- 
gion, the external shear they produce is remarkable and deter- 
mines an expansion of the lens' caustics. 
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Fig. 10. Lensing cross section for arcs with length-to-width ratio larger than 7.5 as a function of the smoothing length. Solid, 
dashed and dotted lines refer to smoothing adopting the asymmetric, the elliptical and the axially symmetric background models, 
respectively. The curves are normalised to the corresponding cross section for smoothing length equal to zero. Shown are the 
results for the cluster projections along the x- (left panel), y- (central panel), and z-axes (right panel). The critical lines in the 
three projections are contained in rectangles sized (395 x 105) h^\pc^, (455 x 205) h^\pc^ and (262 x 87) h^\pc^, respectively. 



The lensing cross sections for large and thin arcs are 
computed as described in detail in several previous papers 
dMeneghetti et al.ll2005l see e.g). We consider here the cross 
sections for arcs whose length-to-width ratio exceeds a thresh- 
old (L/W)min = 7.5, and refer to these arcs as giant arcs. 



f/'i 
ill 



I = 1 1 74 h-' kpc 

I = 704 h-i kpc 

I = 352 h-i kpc 

I = 235 h-i kpc 



Fig. 9. Lens' caustics obtained after removing substructures 
from the region outside a given radius, as given by the labels 
in the figure. The side length is 500 /z 'kpc comoving, corre- 
sponding to ~ 1'. The figure refers to the cluster projection 
along the x axis. 

6.2. Lensing cross sections 

Aiming at quantifying the differences between the strong- 
lensing efficiency of clusters with different amounts of ellip- 
ticity, asymmetries and substructures, we focus on the statisti- 
cal distributions of the arc length-to-width ratios. Indeed, the 
distortion of the images of background galaxies lensed by fore- 
ground clusters is commonly expressed in terms of these ratios. 

The efliciency of a galaxy cluster for producing arcs with a 
given property can be quantified by means of its lensing cross 
section. This is the area on the source plane where a source 
must be placed in order to be imaged as an arc with that prop- 
erty. 



6.2.1 . A particular case: the cluster gShr 

The impact of ellipticity, asymmetry and substructures on the 
lensing cross section for giant arcs obviously depends on the 
particular projected mass distribution of the lens. Large diff'er- 
ences can be found even between different projections of the 
same cluster. As an example, we show in Fig. ^| the lens- 
ing cross sections for giant arcs as a function of the smooth- 
ing length for the three projections of cluster gShr. Results are 
shown for all the smoothing methods described earlier. The 
cross sections have been normalised to that of the unsmoothed 
lens. The horizontal lines in each plot indicate the limiting 
values reached when the surface density maps are completely 
smoothed. Three diff'erent realizations of background source 
distributions were used to calculate the errorbars. 

As expected, the lensing cross sections decrease as the 
smoothing scale increases. The decrement is generally rapid 
for small smoothing lengths, then becomes shallower. 

The differences between the three projections are large. 
When smoothing with an elliptical background model, max- 
imal variations of the cross section of the order of 40% are 
found for the projections along the x- and y-axes. For these 
two projections, circularising the surface-mass distributions re- 
duces the cross section by ~ 85 - 90%. However, while for the 
projection along the .ic-axis smoothing using the asymmetric 
background model reduces the lensing cross section by 35%, 
for the projection along the y-axis the cross section becomes 
only ~ 20% smaller. The dififerences between these two pro- 
jections can be explained as follows. First, as discussed ear- 
lier in the paper, when projected along the x-axis, the clus- 
ter has important substructures close to its centre. This is ev- 
ident in Fig. Is] substructures are significant at radii between 
~ 400 - 800 /i 'kpc. On the other hand, when projected along 
the y-axis the cluster exhibits significant substructures only at 
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g8: low resolution vs high resolution average over four low-resolution clusters 




' 1 ' ' 1 1 

0.1 0.2 0.3 0.4 0.5 

Smoothing length [Mpo/h] 

Fig. 11. Comparison between the low and the high resolution 
version of cluster g8. The lensing cross section for arcs with 
length-to-width ratio larger than 7.5 averaged over three or- 
thogonal projections of the same cluster are plotted versus the 
smoothing function. 

larger radii, > 1 /; 'Mpc. Since strong lensing occurs in the 
very inner region of the cluster, the impact of substructures 
close to the centre is larger than that of substructures farther 
away. Second, in the projection along the y-axis, the twist of the 
iso-density contours and the variations of their ellipticity are 
significantly larger than for the projection along the x-axis (see 
FigS- Therefore, while for the projection along the x-axis the 
deviation from a purely elliptical mass distribution is mostly 
due to the effects of substructures, in the projection along the 
y-axis it is due to both substructures (~ 20%) and asymme- 
tries (~ 25%). Asymmetries which are due to the presence 
of large-scale density fluctuations distort the isodensity con- 
tours which are elongated in some particular direction, vary- 
ing their ellipticity and position angle. Such large-scale modes 
contribute to the shear field of the cluster, pushing the critical 
lines towards regions of lower surface density and increasing 
their size. Consequently, the strong lensing cross section also 
increases. 

As shown in Sect. |5] when projected along the z-axis, the 
cluster appears rounder Consequently, a smooth axially sym- 
metric representation of this lens which conserves its surface 
density profile has a lensing cross section for giant arcs which 
is only 50% smaller than that of the original cluster. Smoothing 
using the asymmetric or the elliptical background models is 
equivalent and leads to a reduction of the lensing cross section 
by ~ 30%. The absence of significant differences between these 
two smoothing schemes indicates that asymmetries play little 
role in this projection, while the large substructure observed at 
~ 8OO/1 'kpc from the centre has a significant impact on the 
lensing properties of this lens, even being at relatively large 
distance from the region where strong lensing occurs. 

The smoothing length for which the curves converge to the 
values for the completely smoothed maps tell us the character- 
istic scale of cluster substructures which is important for lens- 
ing. In those projections where localised substructures play an 
important role, i.e. in the projections along the x- and the z- 
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Fig. 12. Mean lensing cross section for arcs with length-to- 
width ratio larger than 7.5 of four low-resolution clusters as 
a function of the smoothing length. Solid, dashed and dotted 
lines refer to smoothing adopting the asymmetric, the ellipti- 
cal and the axially symmetric background models, respectively. 
The critical regions of the lenses have maximal radii in the 
range ~ 100 - 250/j 'kpc. 

axes, the relevant scales are smaller (< 100 - 300/;"'kpc), 
while for the projection where asymmetries are more relevant 
they are larger (< 400/i"'kpc). Converting these spatial scales 
into the corresponding mass scales in not an easy task, espe- 
cially because we are dealing with substructures in two di- 
mensions. Tentatively, we can assume that the substructures 
are spherical and their mean density corresponds to the virial 
overdensity Ay(z). For z = 0.3, in the cosmological framework 
where our simulations are carried out. A,, ~ 123. Then, the 
above mentioned spatial scales correspond to masses between 
~ 4 X 10" and ~ 2 X 10' V Mq. 

The three cluster projections whose lensing properties were 
discussed above were carried out along the three orthogonal 
axes of the simulation box. In general, these axes do not co- 
incide with the cluster's principal axes because it is randomly 
oriented with respect to the simulation box. Thus, the roundest 
and the most elliptical cluster projections that we have studied 
are not necessarily the roundest and the most elliptical possi- 
ble projections, respectively. In the case of g^hr, however, the 
principal axes do not differ substantially from those of the sim- 
ulation box. The cluster turns out to be prolate with axis ratios 
/1//2 ~ 1.9 and h/h ~ I I- When projected along the major 
principal axis, i.e. in its roundest projection, the ellipticity in 
the central region is slightly smaller than in the projection along 
the z-axis, varying between 0.1 and 0.2. When projected along 
the two other principal axes, the cluster has ellipticity and twist 
profiles very similar to those for the projections along the x- 
and the y-axes. For these reasons, the differences between the 
strong lensing cross sections of the purely elliptical and of the 
axially-symmetric smoothed models are modest in the round- 
est projection, even smaller than for the previously discussed 
projection along the z-axis. Indeed, we find that the ellipticity 
accounts for only 10% of the lensing cross section in this case. 
When projected along the other two principal axes, the impact 
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Fig. 13. Lensing cross section for arcs with length-to-width ra- 
tio larger than 7.5 as a function of the minimal equivalent radius 
containing substructures. Results are shown for all three clus- 
ter projections. The curves are normalised to the cross section 
of the cluster containing all its substructures, corresponding to 
/ - CO. The critical regions of the lenses have maximal radii in 
the range ~ 100 - 250/j-ikpc. 

of the ellipticity is similar to that for the projections along the 
X- and y-axes. 

6.2.2. Mean lensing cross sections 

The example of cluster ^8hr shows that, depending on the par- 
ticular configuration of a lens, the impact of ellipticity, asym- 
metry and substructures can be substantially different in differ- 
ent clusters. Nevertheless, we can try to estimate what is the 
statistical impact of all these factors. For doing that, we repeat 
the analysis shown for the cluster §8hi on our sample of clus- 
ters simulated with lower mass resolution. Among them, we 
analyse the lensing properties of the low-resolution analogue 
of cluster §8hi . When we compare the sensitivity to smooth- 
ing of the low- and the high-resolution versions of the same 
cluster, we do not find significant differences between them. 
Fig. ^2 shows how the lensing cross section for giant arcs 
changes as a function of the smoothing length for all three 
smoothing schemes applied. Each curve is the average over the 
three independent projections of the clusters. The thick and thin 
lines refer to the high- and low-resolution runs, respectively. 
Considering that, as discussed in Sect. |2 the two simulations 
have quite different mass distributions, the differences shown 
here, which are still within the eiTor bars, have little signifi- 
cance. This suggests that our results are not affected by prob- 
lems of mass resolution of the numerical simulations. 

In Fig.^lwe show the variation of the lensing cross section 
vs. smoothing length averaged over all the low-resolution clus- 
ters which we have analysed. For each cluster, we use the three 
projections along the x, y and z axes, i.e. 12 lens planes in total. 
The results are shown for all three smoothing schemes adopted. 
The curves show that on average removing the substructures 
from the clusters reduces their lensing cross section by about 
30 ~ 35%. Removing asymmetries, i.e. transforming the clus- 



ter mass distributions to purely elliptical, further reduces the 
lensing cross section for giant arcs by ~ 10%. If also ellipticity 
is removed, the mean lensing cross section becomes ~ 20% of 
that of the unsmoothed lenses. The typical scales for the sub- 
structures which mostly affect the lensing properties of their 
host halos are < 150/z"'kpc (M < 10'2/z"' Mo). Note that this 
does not mean that larger substructures do not affect the lensing 
cross sections: simply, they are less abundant. We verified that 
only one of the clusters in our sample (§72) is undergoing a ma- 
jor merger with a massive substructure {M^ub ~ 4x10'"*/!"' Mq) 
at z = 0.3. Note also that the largest scale sub-halos con- 
tribute also to the asymmetry of the projected mass distribu- 
tions. This means that smoothing further using the asymmetric 
background model does not remove these large substructures 
completely. When smoothing using the background elliptical 
and axially-symmetric models the smoothing length at which 
the lensing cross section approximates that of the completely 
smoothed lenses is slightly larger, because larger-scale contri- 
butions to the surface density fields must be removed. 

6.2.3. Location of important substructures 

We now investigate what is the typical location of substruc- 
tures which are important for lensing. By ray- tracing through 
the mass distributions obtained after removing substructures 
from outside a given equivalent radius, as discussed at the end 
of Sect|4] we find that the strong lensing efficiency of clusters 
is sensitive to substructures located within a quite large region 
around the cluster centre. For demonstrating this, we show in 
Fig.ll3lhow the lensing cross section changes as a function of 
the minimal radius containing substructures. The cross sections 
are again normalised to those of the unsmoothed lenses. 

In the projection along the y-axis of cluster §8hi (short 
dashed line), we note that the lensing cross section decreases 
quickly when removing substructures outside an equivalent ra- 
dius of ~ 1 /i 'Mpc. The lensing cross section for giant arcs 
is already reduced by ~ 10% when the minimal equivalent ra- 
dius containing substructures is ~ 800/j^'kpc. In fact, in this 
projection there are two large substructures at distances be- 
tween 1 and 1.2/i"'Mpc from the cluster centre, which seem 
to influence the strong lensing efficiency of this lens. Note that 
the critical lines in this projection of the cluster extend up to 
~ 200/1 'kpc from the cluster centre. In the other two projec- 
tions of the same cluster (long dashed and dotted lines), where 
large substructures are located closer to the centre, a similar 
decrement of the lensing cross section is observed at much 
smaller equivalent radii, between 300 and 450 /z^'kpc. 

When averaging over the low-resolution cluster sample, we 
still find that the lensing cross sections start to decrease when 
substructures outside a region of equivalent radius ~ 1 /i^'Mpc 
are removed from the clusters. While the minimal radius con- 
taining substructures is further reduced, the cross sections con- 
tinue to become smaller. The evolution is initially shallow. A 
reduction of ~ 15% is observed at a minimal equivalent ra- 
dius ~ 300/i^'kpc. If substructures are removed from an even 
smaller region around the centre of the clusters, the decrement 
of the lensing cross sections becomes faster. The critical re- 
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Fig. 14. Effects of substructures on scales < 50/z^'kpc on the 
morphology of gravitational arcs. In the left panels shown are 
the arcs formed out of two different sources lensed by the un- 
smoothed mass distribution displayed in the top left panel of 
FigEl In the right panels, shown are the corresponding images 
obtained by using as lens the mass distribution given in the top 
right panel of Fig|2l which has been obtained by smoothing 
with a smoothing length of 47 /z 'kpc. In all panels, the critical 
lines are drawn for comparison. The top and the bottom panels 
have sizes of 38" x 57" and 36" x 56", respectively. 



gions of the lenses in our sample have maximal radii in the 
range - 100 - 250/!"'kpc. 

This shows that substructures close to the cluster centre are 
the most relevant for strong lensing, but substructures located 
far away from the cluster critical region for lensing also have a 
significant impact on the cluster lensing cross sections. 

6.3. Arc shapes, locations and fluxes 

Small changes in the positions of the caustics and therefore in 
the positions of the critical lines can have huge consequences 
on the appearance and location of gravitational arcs. In order 
to describe these eff'ects we compare here the characteristics 
of the images of the same population of background sources 



Fig. 15. Example of gravitational arc shifted by substructures. 
The size of the each frame is 27" x 34". Left panel: simulation 
including substructures. Right panel: simulation performed af- 
ter smoothing with a smoothing length of 47 A 'kpc. 



lensed both with the unsmoothed projected mass distributions 
of the numerical clusters and with weakly smoothed versions 
of them. We smooth using the asymmetric model and using a 
smoothing length of 47 /i^'kpc (M ~ 5 x 10"Vi"' Mq), not sig- 
nificantly exceeding the scale of galaxies in clusters. The aim 
of this discussion is to show that even relatively small substruc- 
tures may play a crucial role in determining the appearance of 
gravitational arcs. 

Sources are first distributed around the caustics of the un- 
smoothed lens following the method discussed in Sect.|3l Then, 
the same sources are used when ray-tracing through the sur- 
face density maps from which substructures are removed. Each 
source conserves its position, luminosity, ellipticity and orien- 
tation, allowing to directly measure the effects that removing 
substructures and asymmetries has on several properties of the 
same arcs. For this experiment, we use an extended version of 
our ray-tracing code which includes several observational ef- 
fects, like sky brightness and photon noise, allowing to mimic 
observations in several photometric bands. We assume that the 
sources have exponential luminosity profiles and shine with a 
luminosity in the B-band Lb - 10'" L©. We simulate exposures 
of 3ksec with a telescope with diameter of 8.2m (VLT-like). 
The throughput of the telescope has been assumed to be 0.25. 
The surface brightness of the sky in the B-band has been fixed 
at 22.7mag per square arcsec. In this ideal situation, no see- 
ing is simulated. The effects of all these observational effects 
on the morphological properties of arcs will be discussed in a 

forthcoming paper. 

In earlier studies. iMeneghetti et alJ ll200rih and lFlores et alJ 

(|2000) showed that the impact of galaxy-sized cluster subha- 
los on the statistical properties of gravitational arcs with large 
length-to-width ratios is very modest. These results are con- 
firmed in the present study. As shown in the previous section, 
the lensing cross sections for long and thin arcs decrease by 
~ 20% when smoothing the cluster surface densities on scales 
~ 50/1 'kpc. On smaller scales the decrement is only of a 
few percent. However, the morphology of individual arcs is 
strongly affected in several cases. Arcs can become longer or 
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shorter, thinner or thicker. In other cases, more dramatic mor- 
phological changes are found. For example, cluster galaxies lo- 
cally perturb the cluster potential such as to break long arcs, 
while in other cases the opposite effect occurs. Two examples 
are shown in Fig^^ In the upper panels, the same source is 
imaged as two short arcs or as a single long arc when the 
unsmoothed (left panel) and the smoothed lens (right panel) 
are used, respectively. The lenses displayed in the top panels 
of Fig. 12 have been used for these simulations. In the bottom 
panels, a single long arc becomes a smaller arclet in absence 
of substructures. Super-imposed on each graph are the critical 
lines. They tend to wiggle around individual substructures in 
the left panels, while they are more regular in the right panels. 
Substructures slightly shift the high-magnification regions of 
a cluster relative to the background sources, inducing remark- 
able changes in the shape of their images and in their multiplic- 
ity. For the cluster projection used in this example, the image 
multiplicity is increased for ~ 21% of the sources producing 
arcs longer than 5", when smoothing is applied, indicating that 
long arcs break up. On the other hand, for ~ 10% of them the 
image multiplicity decreases, showing that the caustics shrink 
and sources move outside of them. Consequently the number 
of their images decreases. 

In several cases, substructures are also responsible for sig- 
nificant shifts of the positions of gravitational arcs. An exam- 
ple is shown in Fig. [21 The size of each frame is 27" x 34". 
The morphological properties of the arc in the two simulations 
are almost identical. The arc length is ~ 11", the arc width 
is ~ 0.6". The luminosity peak of the arc, which we use for 
measuring the shift, is moved towards the bottom left corner 
of the frame by ~ 8.5", when substructures on scales smaller 
than 47 h \pc are smoothed away. Similar cases are frequent. 
For the lens used in this example, ~ 27% of the long arcs 
(length > 5") found in the simulation including substructures 
are shifted by more than 5" after smoothing. About ~ 4% of 
them is shifted by more than 10". From this analysis, long arcs 
which split into smaller arclets are excluded. 

Finally, substructures affect the fluxes received from the 
lensed sources. The histogram in Fig. ^] (solid line) shows 
the probability distribution function of the differences AB = 
B-Bsm between magnitudes of arcs with length > 5" measured 
in the simulations where the unsmoothed and smoothed lens 
projected mass maps were used as lens planes, respectively. 
The analysis is restricted to arcs whose length exceeds 5" in the 
simulation containing all substructures. Some arcs are magni- 
fied, some others demagnified by the substructures. The maxi- 
mal variations in luminosity correspond to AB 2.3 H- -1-2.4. 

The distribution is slightly skewed towards the negative val- 
ues, indicating that in absence of substructures arcs tend to be 
less luminous. In fact, substructures contribute to magnify the 
sources, as discussed in Sect. 16. II 

Sources of different size are expected to be differently sus- 
ceptible to the substructures. The dashed and the dotted lines 
in Fig. ^] respectively show how the probability distribution 
function of AB changes when the source size is increased or de- 
creased by a factor of two compared to the original source size 
used in the simulations. As expected, larger sources are less 
sensitive to perturbations by small substructures in the lenses. 
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Fig. 16. Probability distribution function of the differences be- 
tween arc magnitudes in simulation including and excluding 
substructures on scales < 47/!"'kpc. Results are shown for 
three different source sizes. See text for more details. 



Similar results were found for some other cluster models. 
For other lenses, the impact of the substructures on the proper- 
ties of individual arcs is even stronger 

The observed arc shifts have tangential and radial compo- 
nents. Generally, the tangential shifts are larger than the ra- 
dial shifts. However, when large substructures located close 
to the critical regions of clusters are smoothed away, signif- 
icant radial shifts are possible, given that the relative size of 
the critical lines changes dramatically. In Fig. [n] the radial 
shifts of long arcs (length > 5") is plotted versus the tangen- 
tial shifts. Different symbols are used to identify arcs produced 
by different numerical clusters. As anticipated, for the major- 
ity of the arcs produced by the clusters ^1, g8, g^hr and ^51 
the radial shifts are within few arcseconds, while tangential 
shifts of 10" and more are frequent. On the other hand, the 
arcs produced by the cluster gl2 have significantly larger radial 
shifts. As mentioned above, g72 is undergoing a major merger 
and a secondary lump of matter occurs near the cluster centre. 
The cluster critical lines, along which arcs form, are elongated 
towards it. When moderate smoothing is applied, the impact 
of the merging substructure is attenuated and the critical line 
shrinks substantially. Thus, the arcs move towards the centre 
of the cluster and their morphology and flux are also strongly 
affected. 

Some results for all the cluster models we analysed are 
summarised in Tab^ AH of these effects might have an enor- 
mous impact in lensing analysis of clusters, in particular when 
modelling a lens by fitting gravitational arcs. These results 
show that any substructure on scales comparable to those of 
galaxies should be included in the model in order to avoid sys- 
tematic errors. This problem will be addressed in detail in a 
following paper, in particular regarding the possible biases in 
strong lensing mass determinations. However, by making the 
wrong assumption of axial symmetry, we can approximately 
estimate the errors due to the radial shifts of the arcs. For axi- 



14 



M. Meneghetti et al.: Arc sensitivity to cluster ellipticity, asymmetries and substructures 




Fig. 17. Distribution of long arcs (length > 5") in the plane ra- 
dial (AR) vs. tangential {RA<p) shift. Different symbols identify 
arcs produced by different clusters. The arcs produced by the 
cluster g72, which is experiencing a major merger, are given by 
the small filled squares. 

ally symmetric lenses, the mean convergence within the critical 
line is ^(< Xc) = 1. The mass within Xc is then 



where 

c2 D, 



4;rG DiDu 



(18) 



(19) 



is the critical surface mass density and Xc is in physical units. 
We assume that the position of an arc traces the position of the 
critical line. Then, if an arc distance from the centre changes 
from R to R', the relative variation of the mass inferred from 
strong lensing is 



AM _ M-M' _R^ 



A' 



R'- 



(20) 



The distribution of such AM/M, as derived from the radial 
shifts displayed in Fig[^ is shown in Fig. ^] Without suit- 
ably modelling the effects of substructures, the typical errors in 
mass determinations are within a factor of two, but larger er- 
rors are also possible. Since substructures generally contribute 
to enlarge the critical lines, a larger mass within the critical 
line would be required in order to have an arc at the observed 
distance from the cluster centre. 

Note that even substructures far away from the cluster cen- 
tre are important. For example, keeping the inner structure of 
the projection along the y-axis of cluster ^8hi unchanged, while 
removing the big substructures at distances > 1 /z^'Mpc, we 
find that more than ~ 50% of the long arcs are shifted by at 
least 5". Moreover, image multiplicity increases for 26% and 
decreases for 8% of the sources respectively. 

If relatively small substructures can alter many of the prop- 
erties of gravitational arcs, even asymmetries may be relevant. 
As noted earlier the projection along the y-axis is the most 
asymmetric of g^hr- Comparing the properties of arcs lensed 
by the smoothed asymmetric and elliptical models of this lens, 
we find significant shifts in the location of about 45% of the 
long arcs. For ~ 20% of the sources producing long arcs, the 
multipUcity is changed. 
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Fig. 18. Distribution of the relative variations of mass determi- 
nations from strong lensing, assuming axial symmetry and that 
the arc position trace the location of the critical lines. 

7. Conclusions 

In this paper we have quantified the impact of several prop- 
erties of realistic cluster lenses on their strong lensing ability. 



Table 1. Effects of substructures on gravitational arcs. Column 
1: cluster name; column 2: projection; column 3: percentage 
of sources whose image multiplicity increases; column 4: per- 
centage of sources whose image multiplicity decreases; column 
5: percentage of long arcs (I > 5"), whose positions result to 
be shifted by more than 5" when substructures are smoothed 
away; column 6: maximal variations of magnitudes of long 
arcs. 



Cluster 


proj. 


inc. mult. 


dec. mult. 


shift > 5" 


AS 






[%] 


[%] 


[%] 






X 


21.2 


9.8 


26.6 


-2.3/+2.4 




y 


19.7 


3.9 


11.1 


-1.7/+0.8 




z 


23.7 


10.1 


15.5 


-2.5/+ 1.6 




X 


47.0 


1.0 


24.0 


-1.9/+1.1 




y 


26.1 


0.0 


20.0 


-1.7/+0.7 




z 


24.5 


0.0 


29.5 


-2.5/+0.8 




X 


24.0 


4.0 


27.0 


-1.7/+1.0 




y 


35.4 


9.5 


51.1 


-2.1/+2.1 




z 


31.2 


6.1 


28.2 


-2.5/+0.5 




X 


33.1 


6.2 


44.7 


-2.3/+ 1.7 




y 


35.6 


5.1 


65.7 


-2.5/+ 1.3 




z 


39.1 


6.2 


54.3 


-1.5/+1.3 




X 


36.0 


4.0 


79.4 


-2.5/+2.1 




y 


62.5 


0.0 


57.1 


-2.4/+0.0 




z 


22.5 


0.0 


70.9 


-2.4/+ 1.7 


mean 




32.1 


4.4 


40.3 


-2.1/+1.3 
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In particular, our goal was to separate the effects of substruc- 
tures, asymmetries and ellipticity. For doing that, we analysed 
the lensing properties of one numerical cluster simulated with 
very high mass resolution. In addition, we studied four other 
clusters obtained from N-body simulation with a lower mass 
resolution. 

Each cluster was projected along three independent direc- 
tions. For each projection, we constructed three completely 
smoothed versions. Each of them conserves the mean surface 
density profile of the mass distribution of the cluster. However, 
the first reproduces the variations of the ellipticity and of the 
position angle of the isodensity contours as functions of the 
distance from the centre; the second has elliptical isodensity 
contours with fixed ellipticity and orientation; the third is an 
axially symmetric model. 

The lensing properties of the numerical clusters, of their 
smoothed analogues and of several intermediate versions were 
investigated using standard ray-tracing techniques. 

Our main results can be summarised as follows: 

- Substructures, asymmetries and ellipticity contribute to 
increase the ability of clusters to produce strong lens- 
ing events. Substructured, asymmetric and highly elliptical 
clusters produce more extended high magnification regions 
in the lens plane where long and thin arcs can form. Indeed, 
substructures, asymmetries and ellipticity determine the lo- 
cation and the shape of the lens caustics around which 
sources must be located in order to be strongly lensed by 
the clusters. 

- The impact of substructures, asymmetries and ellipticity on 
the lensing cross section for producing giant arcs is differ- 
ent for different lenses. The lensing properties of the most 
symmetric clusters appear to be particularly influenced by 
the substructures. On the contrary, substructures are less 
important in asymmetric lenses. 

- On average, we quantify that substructures account for ~ 
30% of the total cluster cross section, asymmetries for ~ 
10% and ellipticity for ~ 40%. 

- The substructures that typically contribute to lensing are on 
scales < 150 - 200/i^'kpc. Assuming a virial overdensity 
of ~ 123 for z - 0.3, this corresponds to mass scales of the 
order of ~ lO'^/;"' M©. Substructures on larger scales are 
not as frequent in our cluster sample, but, if present, they 
can boost significantly the lensing cross section (see e.g. 
Ijorri et al. 2004; Meneahetti et al. 2005). 

- Substructures play a more important role when they are lo- 
cated close to the cluster centre. However, the lensing cross 
section for giant arcs is sensitive to substructures within a 
wide region around the cluster core. In particular, our sim- 
ulations show that the sensitivity to substructures far from 
the centre is particularly high in those clusters whose inner 
regions are unperturbed. In these cases, the loss of strong 
lensing efficiency due to removing the substructures from 
the clusters is correlated with substructures within a region 
of ~ 1 /i^'Mpc in radius; on the contrary, clusters contain- 
ing substructures in the inner regions are "screened" against 
external perturbers. 



- Even small substructures (Z < 50kpc, M < 5 x 10"^/i"' Mq) 
influence the appearance and the location of gravitational 
arcs. The perturbations to the projected gravitational po- 
tential of the cluster induced by the substructures alter the 
multiplicity of the images of individual sources. Moreover, 
they change the morphology and the flux of the images 
themselves. Finally, they can shift the position of arcs with 
significant length to width ratios by several arcseconds on 
the sky. 

These results highlight several important aspects of strong 
lensing by clusters. First, any model for cluster lenses cannot 
neglect the effects of asymmetries, ellipticity and substructures. 
Clusters which may appear as relaxed and symmetric, for ex- 
ample in the X-rays, are potentially those which are most sen- 
sitive to the smallest substructures, located even at large dis- 
tances from the inner cluster regions, critical for strong lens- 
ing. Even subhalos on the scales of galaxies can influence the 
strong lensing properties of their hosts and alter the shape and 
the fluxes of gravitational arcs. Therefore, if the lens modelling 
is not caiTied out at a very high level of detail, it may result in 
being totally incorrect. 

Second, the high sensitivity of gravitational arcs to devi- 
ations from regular, smooth and symmetric mass distributions 
suggests that strong gravitational lensing is potentially a power- 
ful tool to measure the level of substructures and asymmetries 
in clusters. Since, as we said, the sensitivity to substructures 
is higher in the case of more symmetric lenses, we conclude 
that dynamically active clusters, like those undergoing major 
merger events, should be quite insensitive to "corrugations" in 
the projected mass distribution but highly sensitive to asymme- 
tries. Arcs could then be used to diagnose mergers in clusters. 
Conversely, substructures should become increasingly impor- 
tant for the arc morphology as clusters relax. Then the level 
of substructures in clusters should be quantified by measuring 
their effect on the arc morphology. This is particularly intrigu- 
ing since measuring the fine structures of gravitational arcs has 
become feasible thanks to the high spatial resolution reached in 
observations from space. 

Third, the strong impact of asymmetries and substructures 
on the lensing properties of clusters and the wide region in the 
cluster where these last can be located in order to produce a 
significant effect further support the picture that mergers might 
have a huge impact on the cluster optical depth for strong lens- 
ing, as suggested in several previous stud ies cTorri et aL..2004: 
iMeneghetti et alJ2005HFedeU et alJ2o'ol . 
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